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Abstract 

There is probably not a black hole in the center of the Sun. Despite this detail, our goal in this work to convince the reader 
that this question is interesting and that work studying stars with central black holes is well motivated. If primordial black 
holes exist then they may exist in sufficiently large numbers to explain the dark matter in the universe. While primordial 
black holes may form at almost any mass, the asteroid-mass window between 10716 — 107! Mọ remains a viable dark 
matter candidate and these black holes could be captured by stars upon formation. Such a star, partially powered by accretion 
luminosity from a microscopic black hole in its core, has been called a ‘Hawking star.’ Stellar evolution of Hawking stars is 
highly nontrivial and requires detailed stellar evolution models, which were developed in our recent work. We present here 
full evolutionary models of solar mass Hawking stars using two accretion schemes: one with a constant radiative efficiency, 


and one that is new in this work that uses an adaptive radiative efficiency to model the effects of photon trapping. 
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1 introduction 


Asking whether the Sun, or any star, has a central black hole 
is equivalently a question about dark matter. While the evi- 
dence for its existence is now overwhelming, the fandamen- 
tal nature of dark matter has eluded explanation. It is clear 
from observations of big bang nucleosynthesis, the cosmic 
microwave background, large scale structure, cluster lens- 
ing, and galaxy rotation curves in spiral galaxies and ve- 
locity dispersions in elliptical galaxies, among a long list 
of other observations, that the majority of the matter in the 
universe simply does not interact electromagnetically (for a 
review see Young 2017). Notably, these observations are in- 
dependent of one another and span cosmic time, probing the 
dark matter abundance from the first minutes of the universe 
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to the present, and are all consistent in finding that 85% of 
the mass in the universe is invisible and thus is dark. 

The composition of this invisible matter is unknown and 
candidates are too numerous to enumerate here, but can be 
broadly grouped based on their masses. Microscopic theo- 
ries of dark matter generally consist of new particles beyond 
the Standard Model (see Sect. 27.5 in Zyla et al. 2020). In 
contrast macroscopic candidates may be composites of new 
particles from an entire dark sector, but many candidates 
only require Standard Model physics such as quark nuggets, 
strangelets, and most notably primordial black holes (Al- 
cock et al. 1993; Witten 1984; Hawking 1971). 

Primordial black holes (PBHs), black holes that may have 
formed in the first instants of the universe, are a popular 
dark matter candidate. Stochastic density fluctuations in the 
early universe may have produced overdense regions that 
collapsed under their own self-gravity, producing large num- 
bers of PBHs (Carr and Ktihnel 2020). While those less 
massive than about 10!5 g could have evaporated due to 
Hawking radiation losses by the present age of the universe, 
those more massive than this limit have barely radiated any 
mass due to the strong scaling of Hawking radiation lumi- 
nosity with mass (L « Me. such that black hole lifetimes 
scale like tay X MB): Observational constraints on PBHs 
with monochromatic masses are now quite tight, and restrict 
PBHs as the entirety of the dark matter to a few mass win- 
dows. Notable for this work is the low mass window be- 
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tween 10716 — 107 !0 Mo, comparable to asteroids (Carr and 
Kiihnel 2022). 

Despite the claimed constraints, PBHs remain a com- 
pelling dark matter candidate for several reasons. First, they 
do not necessarily require new physics beyond the Standard 
Model to form (Carr et al. 2021). Second, they have a simple 
natural production mechanism from the density fluctuations 
in the early universe. And lastly, there is now suggestive ev- 
idence of their existence across large mass ranges (see Carr 
et al. 2023, for a review). 

If low mass PBHs exist and are numerous then it is pos- 
sible that stars could capture them (Lehmann et al. 2022). 
The suggestion that stars may harbor central PBHs, includ- 
ing the Sun, is well founded in the literature. In the seminal 
work, Hawking (1971) estimates that up to 107! Mo of 
PBHs could have been captured by the Sun, and this idea 
has long been entertained in the serious literature. For ex- 
ample, before the discovery of neutrino oscillations it was 
suggested by Clayton et al. (1975) that the two-thirds of the 
Sun’s luminosity may be provided by accretion onto a cen- 
tral black hole, rather that fusion, thus explaining the deficit 
of neutrinos. 

Stars are unlikely to capture PBHs after formation 
(Montero-Camacho et al. 2019). A PBH with finite veloc- 
ity and falling from infinity will be accelerated up above 
the escape velocity and rapidly transit the star. The transit 
timescale is small compared to the damping timescales from 
accretion and dynamical friction, so the PBH is highly un- 
likely to become gravitationally bound to the star or system. 
Capture only becomes probable for PBH transits through 
denser objects such as white dwarfs or neutron stars. Ob- 
servational consequences of a neutron star capturing a PBH 
may include fast radio bursts (Abramowicz et al. 2018), for 
example. Even though capture by main sequence stars is un- 
likely, it is far more likely that a PBH could become bound 
to a star upon formation due to the time dependent nature of 
the gravitational potential of a collapsing cloud. The prob- 
lem of capture during star formation is hard, but past authors 
have considered it and shown that it could be quite likely in 
lower mass halos with slower velocity distributions, espe- 
cially those found in dwarf galaxies or possibly the early 
universe (Esser and Tinyakov 2023). Indeed, arguments are 
now abundant in the literature showing that if asteroid-mass 
PBHs exist in large numbers then one should necessarily ex- 
pect a few to be captured by stars (Capela et al. 2013, 2014; 
Montero-Camacho et al. 2019). 

For many readers, intuition from astrophysics will sug- 
gest that a star that captures a PBH will be short lived and 
look nothing like a star during that life. However, we will 
show that stars with very low mass PBHs could be very long 
lived with many surviving their entire main sequence phase. 
Ultimately the evolution is highly sensitive to the accretion 
physics, which is the subject of the following sections. 
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Much of the theory of accretion of stars with central black 
holes was developed for ‘quasi-stars, star-like objects of 
many thousands of solar masses that may have existed in 
the early universe powered by accretion onto a central black 
hole (Begelman et al. 2008; Ball et al. 2011, 2012). More re- 
cently, these accretion models were adapted for stellar scale 
main sequence stars with asteroid-mass PBHs in their cores 
by Bellinger et al. (2023). To distinguish traditional quasi- 
stars from the scenario where a Sun-like star captures a PBH 
proposed by Hawking (1971), Bellinger et al. (2023) pro- 
posed calling these objects ‘Hawking stars.’ 

This work is a primer on Hawking stars and is both a 
companion and extension of our recent work, Bellinger et al. 
(2023). We begin with an estimate of PBH capture rates in 
Sect. 2. We present a brief description of our implementation 
in the stellar evolution code MESA in Sect. 3, and some basic 
theory of black hole accretion in Hawking stars in Sect. 4. 
In Sect. 4.3 we present stellar evolution models for a range 
of seed masses using the accretion scheme of Bellinger et al. 
(2023). In Sect. 4.4, we present a new accretion model devel- 
oped for this work that accounts for photon trapping, along 
with detailed stellar evolution models. We summarize with 
some open questions in Sect. 5. 


2 Primordial black hole capture rate 


In this section we briefly estimate the capture rates of PBHs 
by stars in the Milky Way using the methods of Esser and 
Tinyakov (2023) and Esser et al. (2023). The mean number 
of asteroid-mass PBHs captured by a star is N = fnv(M) 
where f is the PBH fraction of dark matter, 7 is the ‘merit 
factor’ of the host galaxy (see below), and v(M) is a fac- 
tor of order unity that depends on the mass of the star. The 
factor v(M) is determined numerically from simulations of 
contraction of star forming clouds by Esser and Tinyakov 
(2023), and takes into account the amount of time required 
for the PBH to sink into the star. They show that v(M) is 
largely independent of PBH mass between 10715 Mọ and 
107! Mo, but has some weak dependence above this mass. 
It was only determined for stars up to 0.8 Mo, so future 
work is required for more massive stars. 

The merit factor depends on the dark matter density ppm 
and velocity dispersion o of the halo, and dominates the 
capture probability. If PBHs constitute all the dark mat- 
ter (f = 1) then the merit factor directly gives an order-of- 
magnitude probability that a star-forming cloud will capture 
a PBH. The merit factor is 


=( PDM = a 
™=\ t00GeV/em? ) \ Vie 


and assumes that only PBHs in the low-velocity tail are cap- 
tured. For the Milky Way disk, n ~ 1077, so it is highly un- 
likely that the Sun is a Hawking star even if PBHs comprise 
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all the dark matter. However, it does suggest there could 
be of order 10+ Hawking stars in the Galactic disk. In the 
Gaia Data Release 3 one might expect dozens to hundreds of 
Hawking stars (Gaia Collaboration et al. 2023). In the sec- 
tions that follow, we will show how these can be searched 
for. 

Perhaps more importantly, Esser and Tinyakov (2023) 
shows that most stars in ultra-faint dwarfs like Tucana II 
(n = 0.51) and Triangulum II (7 = 0.95) may have captured 
PBHs upon formation. This has incredible potential for con- 
straining PBHs in the asteroid mass window as dark matter, 
but that is only if the evolution and survival times and of 
Hawking stars are well understood. 

We emphasize that more work is needed on PBH cap- 
tures over a larger range of stellar masses and PBH masses 
in order to put strong constraints on PBH dark matter using 
Hawking stars, but these rough order of magnitude scales 
suffice for motivating this work. 


3 Stellar evolution modeling with central 
black holes 


In the sections that follow we present detailed stellar evolu- 
tion calculations of Sun-like stars with central black holes. 
First, we present a brief discussion of our simulation formal- 
ism. 

The presence of a central black hole has cascading ef- 
fects on stellar evolution, causing the internal structure of 
stars with central black holes to slowly diverge from what 
they otherwise would be in the absence of the black hole. To 
study this, we will need a model for the accretion physics 
and feedback between the black hole and the star. 

The current state of the art is MESA, Modules for Stel- 
lar Evolution, a widely used open-source code for simulat- 
ing stellar evolution in 1D (Paxton et al. 2011; Jermyn et al. 
2023). Despite the equations of stellar structure being cou- 
pled and nonlinear, a few simple changes in boundary con- 
dition along with an accretion scheme is sufficient to model 
stars with central black holes, and was developed recently 
by Bellinger et al. (2023). 

Consider how the equations of stellar structure are im- 
pacted by a central point mass. A central mass changes the 
inner boundary condition for integration from M (r = 0) = 0 
to M(r = 0) = Mpz. The black hole of mass Mpu is there- 
fore treated as point-like relative to the star, and this is not 
such a bad assumption. The Schwarzschild radius of a black 
hole is Rs, = 2GMpy ie, with constants G and c. Stellar 
massed black holes are order kilometers in size, while an 
asteroid mass black hole at the evaporation limit is sub- 
angstrom in size. It should therefore be immediately appar- 
ent that the ability of a low mass PBH to grow by accret- 
ing matter from its surrounding star is severely limited when 
compared to a macroscopic stellar black hole. 
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Even in the absence of accretion, this one change in 
boundary changes gravitational gradients throughout the 
star, and especially in the core. This in turn changes the pres- 
sure and density throughout the star. This has further influ- 
ence on the temperature, and thus nuclear reaction rates. It 
is immediately apparent that even a simple change such as 
a small central point mass results in a nonlinear perturba- 
tion to stellar evolution. In fact, such a change could also be 
used to explore how many other families of non-accreting 
MACHO dark matter could impact stellar evolution if cap- 
tured at formation, and should be studied in detail in future 
work. 

In addition to this boundary condition change, accreting 
black holes can radiate strongly. When also including the 
accretion luminosity, these changes become even more pro- 
nounced. Matter falling into the black hole at a rate Mpg 
must radiate gravitational potential energy at a luminosity 
L, as in a typical astrophysical environment, at a rate 


L= ne (2) 
l—e 

where e is the radiative efficiency such that the mass- 
equivalent energy lost from the system is cMpgu. This is 
generally taken to be approximately 0.1 for a non-rotating 
black hole, but may be as large as 0.30 (or 0.42) for a near 
maximally rotating (or maximally rotating) Kerr black hole 
(Thorne 1974; Begelman 2014). However, it could also be 
near 0 if the photons are trapped and unable to diffuse out of 
the infalling matter sufficiently quickly to escape the black 
hole (Begelman 1979). For now we will use € = 0.08, which 
is approximately the energy released falling from infinity to 
the inner-most stable circular orbit of a non-rotating black 
hole (Bellinger et al. 2023). In practice, the assumption of 
constant radiative efficiency is just a first approximation and 
will be revisited later in this work (see Sect. 4.4). 

One sees from Eq. (2) that an accretion scheme requires 
three pieces of physics: an accretion rate Mpu, an accretion 
luminosity L, and a radiative efficiency €. Expressions for 
any two of these then set the third, and many schemes for 
calculating Mpy and L are present in the literature, result- 
ing in an adaptive € and could be simulated in future work 
(Begelman 1979; Flammang 1982, 1984; Marković 1995). 
In the section that follows we will introduce expressions for 
Lor Mpg and use them to construct two accretion models. 


4 Black hole accretion in the solar core 


We now introduce our accretion models and show some stel- 
lar evolution calculations. Throughout this section, we will 
consider the growth of a fiducial 107! Mọ PBH accreting 
material from the solar core. In Fig. 1 we show how the ac- 
cretion rate and luminosity depend on the black hole mass 
in our first model, which is discussed in detail below. 
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Fig. 1 (Top) luminosity and (bottom) accretion rate of a low mass 
black hole accreting from the solar core, assuming a constant radia- 
tive efficiency of € = 0.08. For low black hole masses the accretion is 
limited by the convective luminosity (solid black) while at higher ac- 
cretion rates the radiation pressure from the accretion luminosity limits 
the accretion to an Eddington rate (dashed black); they intersect at the 
Bondi-Eddington transition (BET), Mggr = 2 x 10720 Mo. The ac- 
tual evolutionary path is shaded. At no point does this model accrete 
at the sound speed limited Bondi rate (blue). This may be a useful toy 
model for a rapidly rotating or magnetized star, and is implemented in 
the stellar evolution calculations shown in Fig. 2. 


4.1 Bondi accretion 


A nominal first estimate for the growth rate would use the 
Bondi accretion rate. For a black hole embedded in an in- 
finite medium, the Bondi radius is Rg = 2G Mgg/ c, with 
Cs the sound speed of the medium. At this distance from the 
black hole, the escape velocity is comparable to the sound 
speed and material from the background begins to freely fall 
into the black hole (Bondi 1952). This region defines the 
Bondi sphere containing infalling mass (87/3) pR? with p 
the density of the medium, and is much less massive than 
the black hole. This is taken in our stellar evolution models 
to be a cavity around the black hole from which it accretes, 
and helps to smooth the numerics at the core (Bellinger et al. 
2023). 
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The rate at which matter is accreted is 


dMpy _ 1pG? Mey o 
dt l ` 


RY 


This differential equation is exactly solvable, 


1 


Mpu(t) = —— ; 
Mpy — (1pG?/2)t 


(4) 


with Mpgu denoting the initial mass of the captured PBH. 
Presently p = 150 g cm7? and cs = 5.5 x 107” cm s7! for 
the solar core. One immediately observes that Mpy(t) di- 
verges to infinity at a finite time tg in the future where 


-1 
tg = 1.2 Myr (>) 
150 g/cm? 


Me o 
E at} 550 km/s 


and if correct would immediately suggest searching for stars 
with central black holes is a fool’s errand. 


4.2 Convection-limited bondi accretion 


For an initial 107!? Mo PBH to grow at the rate given by 
Eq. (3) it would also emit 1076 Lọ at our fiducial radia- 
tive efficiency. This is a luminosity comparable to a very 
low mass star being emitted from the volume of the Bondi 
sphere, which is of order 1 mm. The radiation flux produces 
a steep temperature gradient that drives a convective enve- 
lope, stalling the accretion and significantly depressing the 
accretion rate from this theoretical maximum. 

The convective envelope has a theoretical maximum rate 
at which it can transport luminosity away from the black 
hole. In the theory of quasi-stars this convective maximum 
is 


p 


2 
ar (GMpgu)’ . (6) 


Lpg = l6rn 


where 77 is the convective efficiency (0.1) and IT’; is the first 
adiabatic exponent (5/3). If the accretion luminosity exceeds 
this value then the convective shell will dissipate the energy 
through shocks and quench the accretion (Ball et al. 2011, 
2012). 

By equating Eqs. (2) and (6) one can find a differential 
equation of the same form as Eq. (3). While this equation 
still diverges at a finite time, it does so at 


E p(r) U=ej/e\ ¢ 9 
r= M Gyr (PO) ( 15 (5 


x . (7) 
10-2 i) 550 an) 5/3 
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Fig.2 Top: Kippenhahn diagrams of stellar evolution models of Sun- 
like Hawking stars with varying PBH seed masses. Read vertically at a 
given time one can see the structure of the star from the black hole (dark 
grey), the Bondi sphere (light grey), and the H fusing core (red). The 
Bondi sphere and innermost region of the star is convective (slashed) 
due to the accretion luminosity from the black hole. For the least mas- 


It is now clear that Sun-like Hawking stars, with main se- 
quence lifetimes of 10 Gyr, may live through the entire main 
sequence phase with minimal changes to the appearance of 
the star. It is important to note that this is a consequence of 
the relatively constant core sound speed and density during 
the main sequence. As a star approaches the giant phase and 
the core contracts one should expect the rate to drastically 
increase, necessitating stellar evolution codes to study the 
post-main sequence. 


4.3 Eddington accretion 


We now consider higher mass PBH seeds. While we’ve seen 
that 107!? Mo PBHs can be hiding in Sun-like main se- 
quence stars, a more massive 10~!° Mo PBH would seem 
to destroy the Sun in only 0.34 Gyr from Eq. (7). However, 
the radiation pressure from the accretion luminosity may yet 
again extend the life of the star. In this simple first model, 
the accretion rate is not sound speed limited but rather it is 
radiation pressure limited at the Eddington luminosity, 


Lg =4r Č GMp, (8) 
K 


where « is the opacity of the solar core (x = 1.5 cm? g7!). 
The accretion rate transitions from a Bondi rate to an Ed- 
dington rate approximately when Eqs. (6) and (8) are equal. 


sive seed simulated, a sizeable He core (yellow) accumulates before 
the post-main sequence evolution. Bottom: The luminosity due to fu- 
sion (red) rapidly plummets as the accretion luminosity (black) comes 
to dominate the total luminosity of the star (shaded) when the black 
hole mass is a bit greater than about 107! Mo. Masses (blue) are 
given on the right-hand side axis. See also Bellinger et al. (2023). 


The black hole mass at the Bondi-Eddington transition 
(BET) for the solar core is approximately Mprr © 2 x 
10710 Mo with an accretion luminosity of 1075 Lo. More 
generally, the BET depends on both the metallicity and mass 
of the star as these set x, cs, and po. We emphasize that this 
is only a rough first approximation, and the core gas pres- 
sure can also play an important role in regulating the accre- 
tion rate and should be studied using more detailed accretion 
schemes like those of Flammang (1982, 1984). 

When combining Eq. (8) and (2), we obtain a differential 
equation where Meu œx Mpu, so the black hole grows expo- 
nentially as Mgy(t) = Mpgue’ /TE with a growth timescale 
of 


_ K (1 —e)/e 
tm =0.15 Gyr (54) ( -2 | ©) 


We emphasize that this is not a destruction timescale, but 
an e-folding timescale. The black hole grows by an order 
of magnitude in mass on timescales of tg ln 10 ~ 0.35 Gyr, 
which is long yet again. 

Beginning from near the BET (Mgy = io Mo and 
Lgu = 1075 Lo), it takes approximately 1.4 Gyr to grow to 
1075 Mo. Over the next order of magnitude in mass growth 
the luminosity of the black hole grows from about 107! Lo 
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to be in excess of a solar luminosity, thereby modifying the 
star too drastically for these estimates to be of further use. 

In Fig. 2 we show three examples of the evolution of 
an initially Sun-like Hawking star with a seed PBH in the 
core calculated with MESA. For seed masses less than about 
107! Mo, the appearance of the Sun is largely unchanged 
at present age of the Sun (4.6 Gyr) and the effects of the 
black hole are confined to the inner core. Eventually, the 
black hole grows large enough that the accretion luminosity 
dominates over fusion, causing the star to expand and thus 
quenching fusion. At this time the black hole mass is about 
1075 Mo, or roughly a Uranus mass. We consider this a 
useful marker for delimiting the ‘main sequence’ from the 
“post-main sequence.’ 

Once in its post-main sequence phase the Hawking star 
becomes fully convective and slowly swells over multiple 
Gyr to about 10 Ro and about 10 Lo causing them to ap- 
pear as a sub-subgiant star, joining a relatively sparse pop- 
ulation on the HR diagram. This phase of stellar evolution 
is qualitatively insensitive to the initial seed mass, and lasts 
for multiple Gyr in our simulations. Simulations with lower 
constant radiative efficiencies produce similar results, but 
with shorter post-main sequence phases due to the greater 
accretion rates (omitted for length). 

The transition of a Hawking star from its post-main se- 
quence phase to some other high energy transient is not re- 
solved here, as the MESA models fail when the Bondi radius 
reaches the surface of the star (Mgy ~ 107! Mo). Signif- 
icantly more work is required to bridge the late stages of 
evolution in simulations. 

For the lowest mass seed simulated (107 !2 Mo) we find 
that the Hawking star survives for the entire main sequence 
and even succeeds in forming a He core surrounded by a H- 
burning shell. At the onset of He core contraction the black 
hole accretion rate rapidly increases due to the increased 
core density, sound speed, and conductivity. The resulting 
increase in accretion luminosity causes the inner convec- 
tive region to grow outwards, mix the He core with the H- 
burning shell, and quench fusion. In contrast, those models 
with more massive seeds reach the post-main sequence be- 
fore He core formation, so the Hawking star transitions more 
gradually to its post-main sequence phase with a growth 
timescale closer to what is predicted in Eq. (9). 


4.4 Photon trapping and super-Eddington accretion 


In this section we introduce the second accretion model that 
will be considered in this work. The Eddington accretion 
scheme above assumes that photons are not trapped, mean- 
ing that photons can diffuse out from the infalling matter in 
the Bondi sphere faster than the infalling matter can trans- 
port the photons inward. However, it is possible the photons 
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Fig. 3 (Top) luminosity and (bottom) accretion rate of a low mass 
black hole accreting from the solar core assuming the adaptive radia- 
tive efficiency of Eq. (12) and fiducial values for the solar core. As the 
black hole mass grows it more efficiently transports photons inward, re- 
ducing the luminosity. The black hole goes ‘dark’ as it rapidly evolves 
to higher masses by accreting at the Bondi rate given by Eq. (3). 


are trapped, especially in slowly rotating or weakly magne- 
tized stars. As a consequence, one should expect the radia- 
tive efficiency to decrease, especially for higher mass black 
holes, resulting in rapid accretion at late times. 

In Fig. 3 we show a second accretion model using an 
adaptive e€, discussed in detail below. As a brief introduc- 
tion, this model uses € = 0.08 for low black hole masses 
but € rapidly goes to zero above a mass threshold near the 
BET. To reiterate, there are three variables — L, Mpu, and 
e — two of which must be set. In both accretion models we 
set €, leaving one of either L or Meu to be chosen while 
the other is be determined by Eq. (2). The convective limit 
and Eddington limit are imposed by luminosities, while the 
Bondi limit is imposed by the accretion rate (i.e. the maxi- 
mum rate matter can enter the Bondi sphere). Thus, the con- 
vectively limited luminosity and Eddington luminosity are 
unchanged from Fig. 1, while the associated accretion rates 
asymptotically approach infinity with increasing Mgy due 
to the declining €. Meanwhile, the Bondi accretion rate is 
unchanged as this is a limit imposed by the sound speed of 
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the ambient accreted material. So, in contrast with the pre- 
vious accretion model, we see that as the Meu needed to 
maintain the Eddington luminosity grows exponentially, it 
eventually reaches a maximum at the accretion rate provided 
by the Bondi rate. 

We now derive our model in detail. The mean free path of 
a photon is £ = (kp)! =4.4 x 1074 cm for the solar core. 
Assuming a random walk, such a photon has a diffusion co- 
efficient of D, = c£/3. The diffusion coefficient allows us 
to calculate a typical photon diffusion timescale t,, during 
which a photon at speed c taking steps of length £ explores 
(and escapes from) a sphere of radius R in ty = R?/ 6D,. 
Taking R to be the Bondi radius, the photon trapping time is 


2 =4 
EE Mpuy Cs 
ty =1.1 ps 10 
2x 10-10 Mo) \550km/s 


p K 
x . 
150 g cm7? ) (1.5 cm? g7! 


To determine if photons are trapped, we compare this to 
the gravitational free fall time of the infalling matter that 
the photons are diffusing through. This can be found from 
Kepler’s third law by assuming that a vertical free fall is 
equivalent to an elliptical orbit contracted to maximum ec- 
centricity. Such an orbit has a semi major axis of Rg /2, and 
we take the half-orbit time to be the free fall time, 


or = 0.5 us Msn w es (10) 
ee? PN IX 10-10 Mo J \550km/s) ` 


Observe that the photon diffusion timescale and matter 
free fall timescales are comparable near the BET. Since the 
diffusion time scales like Mou while the free fall time scales 
like Mgu, one should expect photons to diffuse out much 
faster than the free fall time for Mgy S Mper while for 
Mpu =, Mper the Bondi sphere should be trapping photons, 
potentially allowing for rapid Bondi-like accretion that is 
only limited by the supply of infalling matter at the Bondi ra- 
dius. This would simultaneously increase the accretion rates 
above the prediction from Eq. (8) while decreasing the ra- 
diative efficiency e€, as the matter-energy that is liberated as 
photons are ultimately consumed by the black hole. 

To incorporate photon trapping in a simple accretion 
scheme we use an adaptive mass-dependent radiative ef- 
ficiency rather than fixing € = 0.08. For a simple model, 
one can assume photons must explore and escape the Bondi 
sphere of radius Rg in a free fall time tr, i.e. they must 
diffuse outward faster than the infalling matter can advect 
them across the event horizon. As photons explore the space 
they become distributed with respect to the origin, with a 
radial distribution that could be taken to be a Rayleigh dis- 
tribution. For small black hole masses, almost all photons 
escape. At greater black hole masses only those in the tail of 
the distribution far from the origin escape. The fraction that 
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are trapped can be taken to be the Rayleigh cumulative dis- 
tribution function 1 — exp(—R3 / (207)) where ø is a scale 
parameter with o? = 6 Dtp for a 3D diffusive random walk 
such that 


Re _ GMpypKk 


202 TCCs 


(1) 


Thus, the adaptive radiative efficiency might be (at fiducial 
values) 


e =0.08exp | —1.15 Mex (12) 
ee eh SDM J |" 


This adaptive radiative efficiency smoothly transitions from 
near 0.08 at M ~ 1071? Mo to near 0 at M ~ 1078 Mo. 

In Fig. 4 we show two stellar evolution models that depict 
the rapid late time growth of the central black hole with pho- 
ton trapping for initial 107!! Mo and 107!? Mo seeds. For 
numerical reasons both the radiative efficiency and accretion 
luminosity have a floor of 1.1 x 107%. At early times when 
e © 0.08, the model with an adaptive radiative efficiency is 
largely indistinguishable from the first model, and for suf- 
ficiently low seed masses (Mgy < 107!* Mo) the star sur- 
vives for its entire main sequence and accumulates a large 
He core. However, as a consequence of this new accretion 
scheme Sun-like Hawking stars do not have an extended 
post-main sequence evolution. 

The stellar evolution models with an adaptive € begin to 
diverge qualitatively from those with a fixed € near the BET. 
When the star has about a million years left to live it tran- 
sitions to a brief Eddington-limited phase, during which the 
accretion luminosity reaches a maximum at about 1075 Lo 
before precipitously dropping with the onset of trapping. For 
example, the central convective region in the model with a 
107!? Mo seed expands out to 1% of the radius of the star, 
but does not succeed in mixing the H burning shell with the 
He core before trapping causes the convective core to shrink 
inwards. 

The black hole mass at the onset of trapping shows weak 
sensitivity to the seed mass. The model with the 10~!* Mo 
seed has accumulated a large He core, so the higher density 
causes the onset of photon trapping at a lower black hole 
mass than the model with a 107!! Mo seed. Since higher 
mass seeds also reach the trapping condition during the H 
burning phase, the late time evolution is the same as model 
with a 1071! Mo seed. 

In the absence of a large accretion luminosity the star 
does not expand, and the black hole grows by about 10° 
over about a thousand years with almost no changes to the 
star’s external appearance. In the final hours of the simu- 
lation, when the black hole mass approaches 0.01 Mo, the 
core temperatures and densities of both models rise by about 
a half order of magnitude and drive rapid increases the rates 
of CNO and pp burning. The growth in fusion luminosity 
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Mpgy=107-''Mo 


Radius (Ro) 


Luminosity (Lo) 
Mg /Me 


Time Until Destruction (yr) 


Fig. 4 (Top) Kippenhahn diagrams of late stage evolution of a central 
black hole for a Sun-like Hawking star including photon trapping. Col- 
ors and shading is the same as Fig. 1. Early times, prior to the transition 
to an Eddington-like accretion, are largely indistinguishable from the 
models with constant € in Fig. 1. The break when about a thousand 
years remain is associated with the onset of trapping, while the final 


can be seen in both models in the bottom panels of Fig. 4 
and is thought to play a role in the flattening of Mgy in the 
final hours. 

At this point the accretion timescale and simulation 
timestep are comparable to the hydrodynamic timescale of 
the Sun and the simulation stops for numerical reasons. 
Therefore, even though the final cannibalism of the star 
proceeds quickly with photon trapping, there may still be 
interesting transients associated with direct collapse due to 
the rapid onset of late-time burning that should be explored 
in future work. While nuclear burning has been largely ig- 
nored in past work studying direct collapse of stars with 
central black holes (e.g. Marković 1995), this work finds 
that it may be important and merits revisiting with detailed 
stellar evolution models. 


5 Summary 


In this work we have presented some elementary theory of 
black hole accretion and used it to calculate stellar evolution 
models of ‘Hawking stars’ — main sequence stars that may 
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break with a few hours remaining is associated with a rapid increase in 
nuclear burning. The time until destruction is taken to be the time until 
the end of the simulation. (Bottom) The accretion rapidly accelerates 
as the accretion re-enters the Bondi regime, and the mass grows rapidly 
to about 107? Mo before the simulation ends for numerical reasons. 


have captured very low mass PBHs in their cores during for- 
mation. 

The two models presented in this work may be best 
thought of as extremes, showing how a high radiative effi- 
ciency or a low radiative efficiency can drastically influence 
the evolution toward a long-lived sub-subgiant phase or a 
rapid direct collapse. It is not immediately clear how suit- 
able either of these models are for describing the real physics 
of microscopic black hole accretion in a main sequence star 
in the modern universe. One might reasonably guess that a 
star could evolve following one of these two qualitative evo- 
lutionary tracks depending on the stellar mass, metallicity, 
rotation, and magnetic field. In this work, we have only ex- 
plored one fiducial star and a few seed masses using very 
simplified accretion schemes. Future theory work will be 
needed to understand accretion in Hawking stars and imple- 
ment these accretion schemes in stellar evolution codes. For- 
tunately, a large body of theory on microscopic black hole 
accretion already exists and could be implemented in MESA 
in the near future (e.g. Begelman 1978, 1979; Thorne et al. 
1981; Flammang 1982, 1984; Marković 1995,?). 

Observationally speaking, either evolutionary outcome is 
interesting. While long lived Hawking stars of about a so- 
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lar mass may spend billions of years as sub-subgiants or red 
stragglers, they may be distinguishable from traditional stars 
via asteroseismology due to the fully-convective nature of 
the star (Bellinger et al. 2023). Meanwhile, rapid direct col- 
lapse provides a mechanism to produce sub-solar mass black 
holes that are otherwise not predicted by traditional evolu- 
tionary channels. As a central black hole does not signifi- 
cantly modify the early phases of evolution, regardless of 
photon trapping at later times, we suggest that many main 
sequence stars could be harboring central black holes. In- 
deed, this work finds that the Sun could be Hawking star 
with a soft constraint on the PBH seed mass being less than 
about 107!! Mo. 

Many open problems must be addressed to make accu- 
rate models of stellar evolution with central black holes; a 
few are discussed below. In rotating magnetized stars, does 
the infalling gas form a disk and drive a jet (Blandford and 
Payne 1982)? Many stars are rapidly rotating and could have 
Kerr parameters that are near maximally rotating. Does the 
black hole get spun up as it accretes, and how does this effect 
the radiative efficiency (Marković 1995)? While disks and 
jets will not be modeled directly in a 1D code like MESA, 
their effects could be mocked up by appropriately chosen 
accretion rates and luminosities. And what is the convective 
efficiency in the Bondi regime (Jermyn et al. 2022)? What 
is the radiative efficiency? Do arguments used to derive the 
radiative efficiency in the continuum limit apply for black 
holes with sub-angstrom Schwarzschild radii? If not, should 
the efficiency be greater or lesser, and at what black hole 
mass do they transition to the canonical efficiency? How 
much mass is lost in winds during the final stages, and what 
is the final black hole mass? These are just a few questions 
that must be addressed. 

Beyond open problems important for the stellar model- 
ing, there are additional questions that must be answered to 
make the results of stellar evolution calculations astrophys- 
ically useful for constraining PBH dark matter. What is the 
capture rate of PBHs in stars? This will depend strongly on 
the PBH velocity distribution in star forming halos, which 
is itself a function a cosmic time. How does the capture rate 
scale with stellar mass? While this work uses the Sun as a 
test case, other Hawking star masses should be explored in 
future work. PBHs may also form in gravitationally bound 
clusters, which further complicates calculations of their in- 
teractions with starforming clouds (Belotsky et al. 2019; 
Gorton and Green 2022). If a PBH becomes bound to a star 
forming cloud, how long does it take to migrate toward the 
center? Are there X-ray transients associated with stellar de- 
struction that are distinct from other accreting binary sys- 
tems, and if so, what are the lifetimes of this transient? An- 
swering these questions is likely beyond the scope of stellar 
evolution models and simple order of magnitude estimates, 
and requires detailed calculations. 
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